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The explosion of a core collapse supernova is approximated by the breakdown of the steady-state 
solution for accretion onto a proto-neutron-star (PNS). We analytically show that as the neutrino 
luminosity exceeds a critical value L c , the neutrinosphere pressure exceeds the hydrostatic limit even 
for an optimal shock radius R. This yields L c oc M 2 T 2 (with logarithmic corrections) and R oc M/T, 
in agreement with numerical results, where M, T are the PNS mass, neutrino temperature. The 
near-critical flow can be approximated as a ballistic shell on top of an isothermal layer. 



Despite decades of research, the revival of the stalled 
shock wave that forms around the collapsed core of a 
massive star is not yet understood. The consequent su- 
pernova explosion is reproduced in general only in low, 
M < 11M Q mass stars [l[, using sophisticated multidi- 
mensional numerical simulations 0]. 

A simplified method of testing for an explosion is to 
examine the conditions under which a steady-state flow 
can exist between the proto-neutron star (PNS) and the 
shock Although this approximation does not cap- 
ture the full complexity of the process, it does reproduce 
the critical behavior: when the neutrino luminosity L 
exceeds a threshold L c , no steady state solution exists, 
corresponding to an instability leading to an explosion. 
However, the nature of the critical behavior, the depen- 
dence of L c upon the flow parameters and the dimension- 
ality, and the nearly constant value of the shock radius, 
R ~ 200 km, are poorly understood (but see Some 
critical criteria have been suggested previously [||-|6| , but 
have not produced generic correct results for L c and R. 

In this Letter we identify the global origin of the critical 
behavior, and analytically derive L c and R for arbitrary 
parameters and equations of state (EoS). We then study 
the flow structure, and show that it is well approximated 
by a simple, single- or double-component model. 

Steady flow model - A stationary, spherically symmet- 
ric flow of mass accretion rate M onto a central mass M 
is governed by the conservation of 



mass, 

momentum, 
and energy, 



Airr pv = M ; 
dv 1 dP _ GM 
dr p dr r 2 
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The variables v, p, T, P and e are the inward velocity, 
mass density, temperature, pressure, and specific internal 
energy of the flow, respectively. Gravity is treated in the 
Newtonian limit, and self- gravity of the accreting layer is 
neglected; G is Newton's constant. The specific energy 
deposition rate q — qn — qc is due to neutrino heating 



minus neutrino (Urea) cooling, as approximated by @] 



QH = k H Lr~ 



qc 



k c T b 



(4) 



where k H ~ 1.5 x 10 20 {L 52 7frJ ) erg s" 1 g"\ k c ^ 
1.4xl0 20 (T/2 MeV) 6 erg s" 1 g _1 , r = lOOrioo km is the 
radial position, and T = 4J~4 MeV is the neutrino tem- 
perature at the neutrinosphere (subscript v), r — r v . The 
single-species neutrino luminosity, L = 10 52 L52 erg s , 
is approximated as black body radiation, so 



L = Airrl (7/16) aT 4 



(5) 



where a is the Stefan-Boltzmann constant. 

We assume nearly pressureless free fall above R, so the 
upstream (subscript it) conditions are determined by 



vl = aGM/R , 



(6) 



where a is a constant of order unity. The flow boundary 
conditions at R are then dictated by the downstream 
(subscript d) properties of a strong shock, given by the 
Rankine-Hugoniot conditions 



VdPd = V u p u 
Pd = Puvl 

h d = e d + Pd/Pd - (vl - v 2 d )/2 



(7) 



Pdv 2 d ; 



where h = e + P/p is the specific enthalpy, and energy 
losses due to dissociation at the shock are neglected. 

The equation system is closed by specifying an EoS for 
the matter and radiation, and by imposing the optical 
depth of the accretion layer Q , 



up dr 
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(8) 



where the specific opacity is given by p| 0] 

k — 9.3 x IQ-^Tl cm 2 g" 1 . (9) 
Alternatively, one may fix the neutrinosphere density 

Pv = p(r v ) = p* (10) 



2 



Q, where p* ~ t/(ktv) ~ 10 11 g cm 3 . This condition 
is less physical but more tractable than Eq. ([5]). 

Solutions to Eqs. (HHTJ) supplemented by either Eq. © 
or Eq. (|10[) are in general found for a given EoS only for 
L smaller than a critical luminosity L C (M, M, T, a). 

Neutrino sphere region. - Near r„, neutrino heating and 
cooling are approximately balanced, qjj — qc, so [5j 



r(»v) 



1/6 



-1/3 _ 



T 



= _^0.79T. (11) 



In this region, the strong cooling and its sensitive, qc oc 
T 6 temperature dependence, indicate that T is strongly 
regulated, such that approximately T(r > r„) ~ To. For 
p of order 10 11 g cm -3 and T of order a few MeV, non- 
relativistic nucleons dominate the pressure, so P ~ fopTo, 
where 6 = fcp/m p , fcp is the Boltzmann constant, and 
to p is the proton mass. The specific kinetic energy and 
its derivative gv = d(v 2 )/dr are small even near R, and 
are completely negligible near r v . Eq. ([2} thus implies 
that p declines roughly exponentially near r v , 



p(r > r„) ~ p(rv)exp 
~ p(rv)exp 



r„<5 



with a scale height r u 6 oc L/(MT 3 ), where 



bT r w 
GM 



0.1 



74^,60 

Mi. 4 



(12) 
(13) 

(14) 



Mi .4 = (M/1.4M Q ), and rv j60 ee (rv/60 km). The neu- 
trinosphere region r < r„(l + (5) thus makes an important 
contribution to the optical depth of the accreting layer, 

r w Kp(r v )r u 5 m Q.§T%rl fi0 p Vf iiM^l , (15) 

where p^u ee [p(r l/ )/10 11 g cm -3 ]. 

Critical behavior.- Next, we analytically show that r 
[Io| and p y are both ( j) maximized at a certain value of R 
that depends to first order only on M/T; and (ii) mono- 
tonically decreasing functions of L for any fixed R (in 
particular, at the maximum). This pinpoints the critical 
behavior. For L < L c there are typically two solutions R 
(only the smaller value is physical 0) for which r = r*. 
But for L > L c , no solution exists because r < r* for all 
i?. Similar behavior is found for p v and p*. 

As gv is small even behind the shock, Eq. ([2]) yields 



Pv-Pd 



W 



R GMp 



bT pu < 1 - exp 



(16) 



R 



where we used Eq. p^|) to approximate p in the near 
vicinity of r u . Note that due to the r~ 2 factor in the 



integrand, this approximation is even better here than it 
is for computing r, e.g., in Eq. (|T5|) . Hence, 



Pd 

bTr 



exp 



i? 



(17) 



Consider hrst the requirement p v — p« , rather than the 
closely related but more complicated r criterion. Eq. (|17[) 
shows that for a given R, p v admits a maximum at 



R = 



GM 
bT \p 



2UM 1A ZT 1 \ P 1 2 5 km 



(18) 



where Ap ee 2.5Ap,2.5 = —dhiPd/dlnR is the radial 
power-law index of the downstream pressure. Further- 
more, for this (or any fixed) value of R, p v is a monoton- 
ically decreasing function of L, 



dpv r _ x dp v Q 
dL v dr v 



(19) 



These two conclusions confirm the critical behavior dis- 
cussed above: criticality corresponds to the R value that 
maximizes p v , and to the r v value for which this maxi- 
mum equals p*. Finally, combining equations (|5|), (|17[) . 
and (|18|) . yields the (p v = p*) critical luminosity, 



7na 

W 3 



GMT jb 



\ P + \n(bT p*/P d ) 



(20) 



Next, we modify these arguments for the r = t* crite- 
rion in Eq. ([8]). For simplicity, we relate r to p„ by in- 
tegrating the approximate form of p(r) in equation (|13p . 
Plugging p^ into Eq. (fT7| then gives 



PdK s - 

k— — e 

GM 



x (i-*) 



1 



L (^-!) 



(21) 



Here, the second term in the square brackets is expo- 
nentially small. Neglecting it and requiring dr/dR = 
reproduces the shock radius in equation (|18l) . The neutri- 
nosphere radius is then found by solving r = r* , yielding 
the critical luminosity 



^ - 2 4/3 



GMT/b 



2W - 



y<e->-PP d GMK/T 
2bT 



(22) 



where W(z) is the Lambert W (branch k = — 1 of the 
product log) function defined by z = We w . 

Comparison to numerical results. - The above analysis 
agrees well with numerical computations, considering the 
level of approximation. To see this, consider a simplified 
EoS, describing nonrelativistic nucleons and relativistic 
electrons of zero chemical potential and degeneracy, 
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(23) 
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FIG. 1: Critical values for the test case as a function of M , 
computed numerically (symbols with thin lines to guide the 
eye) and analytically (thick curves). Shown are L c (left axis), 
computed both for t* = 2/3 (red, filled circles, solid curves; 
Eq. I22[) and for p» = 10 11 g cm -3 criterion (black, empty 
circles, dashed curves; Eg. I20[). and R computed for t* (blue, 
right axis, triangles and dot-dashed curves; Eq. I18[l . 



with a being the radiation constant. Figure [T] shows 
the corresponding values of L c and R for parameters 
{-Mi. 4, 74, a} all being unity (henceforth: the test case). 

At high L (high M), the model accurately reproduces 
the numerical results, including the difference between 
the t = r* and the p — p* criteria. At lower L and M, the 
model increasingly overestimates the numerical results. 
The error is mostly due to an increasing deviation of p(r) 
near r v from the exponential form of Eq. (|12[) . 

For illustration, we fit the downstream pressure by 



P d ~ 5.5 x 10 27 M°fM°"R- 2 52 a - 55 



erg cm 



(24) 



around the test case parameters and Mi ~ 1, where 
i?ioo = (-R/100 km). This is nearly the result ob- 
tained for a fixed shock compression ratio £ ~ 5: Pd = 
(47r)- 1 (aGM) 1 /2(i _ ^)MR- 5 / 2 . It yields 



L 



7.0MlfM [ ^ 22 T 4 2 



erg s 



and 



L 



7AMli°M\ 



0.30-t-: 
'4 



' 2 77 a 017 erg s" 



(25) 



(26) 



These results overestimate their numerical counterparts, 



L^> ~ 5.42 and L)^'' ~ 5.21, by - 30%. Eq. (p]l 
overestimates the numerical i? ~ 225 km by < 10%. 

Isothermal model. - Next we seek the simplest descrip- 
tion of the flow that captures its essential properties and 
critical behavior. Such a model would provide a useful 
tool for analytically testing various generalizations and 
additional physical processes. 

The large contribution of the near-r^ region to r, and 
the tightly controlled T profile in this dense region, moti- 
vate an isothermal model for the flow. For T{r) — const., 



(p.) 



Eq. ([2]) with a negligible kinetic term yields the exponen- 
tial density profile of Eq. §F2§, and Eq. © yields T{r) = 
To to lowest order in S. In a simple isothermal model, 
this behavior is assumed to persist for all r v < r < R. 
Requiring that p v ~ Pd(R) exp[<5 -1 (l — r v /R)] is maxi- 
mized at some R then reproduces Eqs. (TT5|) and (|20l) . but 
with pressure replaced by density, 



R : 



GM 
bT X p 



lira 

24/3 



GMT lb 



X p + ln(p*/pd) 



(27) 



where A p = 



d\n pd/dlnR. Thus, this simple model 
somewhat overestimates the true values of R and L c , but 
approximately reproduces their correct functional depen- 
dence. Similar results are obtained if r, rather than p„, 
is held fixed. 

Isothermal/ballistic profile.- Although the isothermal 
toy model qualitatively reproduces the correct critical be- 
havior, it yields unrealistic flow profiles. Next we general- 
ize the model by introducing an additional, ballistic shell 
lying above the isothermal component. 

The accretion layer has long been recognized as nearly 
hydrostatic, in the sense that the specific kinetic energy 
term gy in the momentum Eq. ([5]) is small. We find 
that in the outer part of the accretion layer, near critical 
conditions, this term is in fact nearly constant, |9 r i>| <C 
v/r; see Figure [5J Such ballistic behavior corresponds, 



by Eq. Q, to p 



and appears in general as L 



approaches L c . Here, heating dominates over cooling in 
the outer part of the flow, so q ~ qn — kuLr^ 2 . 

It is difficult to analytically solve Eqs. ((J)-© even 
in this limit {i.e., neglecting cooling and kinetic terms). 
However, these equations do admit a simple p en r~ 2 
solution if one considers only the relativistic particles, 



Jrel) _ 



GMM 
\2nk H L 



(GM) 2 M 



(28) 



Mathematically, the equations also admit a simple p oc 
r~ 2 solution if one considers only nucleons, albeit with 
a different T = GM/ (3br) profile, and a nonphysical 
/2 < 0. It is therefore not surprising 



P 



(nuc) 



-p 



(rel) 



that the full solution approaches a p oc r behavior. 

Throughout the ballistic region, qu oc r~ 2 is smaller 
than the rate of change of specific gravitational energy, 
GMv/r 2 . The situation is reversed in the isothermal 
region, inward of r g defined by 



Pa = P(r g ) 



GMM 
Airk H L 



(29) 



It is therefore natural to adopt r g as the interface radius 
between isothermal and ballistic shells. The model thus 
becomes 



p(r) ~ \ Pu 6XP L*V 



1)] if 



r < r a 



if r > r a 



(30) 
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FIG. 2: Critical structure of the flow for the test case parame- 
ters with Mi = 1 and L52 = 5.42. Shown are the density (thin 
red curves; left axis) and temperature (thick blue curves; right 
axis) radial profiles according to the numerical solution (solid 
curves) and to the two-component model (dashed curves) of 
Eqs. (f29)l " (p2)) . Also shown are the velocity 11 = 10 9 U9 cm s -1 
(dot-dashed; left axis) and optical depth (upper axis) profiles. 



and 



if r < r, 



T (r/r g )- 3 / 4 iir>r g 



B ' 



T(r) 



where we adopt the T(r) scaling of Eq. ([28 
continuity implies 

GM 



(31) 



and p(r) 



2bTnW 



-l/i GMp„k H L 
Mb 2 T? 



(32) 



The simple model in Eqs. (|29|) - (|32j) is typically accurate 
to within 25% in p, and 15% in T; see Figure [5] The accu- 
racy drops to within a factor of ~ 2 at extreme accretion 
rates (M x < 0.1 or M x > 10). 

Summary. - We identify the global breakdown of steady 
state accretion solutions onto a PNS with the upper limit 



on the optical depth (or p v ) of the accretion layer. This 
limit arises because for L > L c , the neutrinosphere pres- 
sure P v exceeds the hydrostatic limit Pd + W (see Eq. [T6|) 
even for the optimal value of R, which sets the critical 
shock radius. This provides approximate critical crite- 
ria (Eqs. [THl [20l |22|) which agree with numerical com- 
putations, in particular at high M. The accretion layer 
is approximately described as an isothermal component 
of temperature T (Eqs. [TTl - fT4")) , which alone reproduces 
critical behavior, lying underneath a ballistic shell in 
which p oc r~ 2 (Eqs. 



The analysis is quite robust, and the main results are 
independent of the exact forms of the cooling function 
and the EoS. For simplicity, we neglected self-gravity 
and shock dissociation, and assumed Newtonian gravity, 
near-pressureless free fall above the shock, and a constant 
L(r). These assumptions can be relaxed, and the analysis 
generalized, by modifying Pd in Eq. (|17[) or introducing 
corrections to the p(r ~ r v ) profile in Eqs. (fTT]) -(fT4" |) . 

Our results provide useful insights to various aspects of 
the critical behavior. For example, the maximal nature 
of t or p v immediately implies that L c should decrease at 
higher dimensionality, with a higher degree of substruc- 
ture, and in the presence of pulsational phenomena, in 
concordance with the results of 0, H [p^Wla ] . The reason 
is that spatial variations in the Sow, which are inevitable 
at high dimension and are enhanced by substructure and 
pulsation, lead to local variations in the effective values of 
parameters such as L, T, M and a. The larger the vari- 
ations, the more likely it is that some regions in the flow 
would enter the phase space in which the only solutions 
have insufficiently small r or p vi leading to instability. 
Smaller values of L are needed to avoid the onset of such 
local instabilities, diminishing L c . 
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